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ABSTRACT 

Radiative shock waves may be subject to a global thermal instability in which the 
cooling layer and shock front undergo growing resonant oscillations. For strong hy- 
drodynamic shocks, the presence of the overstability depends on the temperature and 
density indices of power-law cooling functions and the specific heat ratio, a, /3 and 
7, respectively. Here, we investigate the stabilising influence of a transverse magnetic 
field by introducing the shock Alfven number, M a as a fourth parameter. We thus 
investigate the stability criteria for both molecular and atomic shocks under a wide 
range of conditions. In particular, we find that all molecular shocks in which the cool- 
ing increases with the temperature (a > 0) are stabilised to the first four modes if 
M a < 20 (fi = 2). For a = —0.5, the first overtone remains stable only for M a < 8. We 
conclude that molecular shocks in the interstellar medium are probably stabilised by 
a transverse magnetic field unless exceptional circumstances arise in which the cooling 
strongly increases as the gas cools. 
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1 INTRODUCTION 

To interpret observations of an astrophysical shock we need 
to know if it can be compared to steady state models or if it 
is subject to an instability. In fact, the cooling immediately 
downstream of a shock front may lead to an overstability of 
the entire radiative shock wave (?). For hydrodynamic flows, 
a linear analysis yields conditions for growing oscillations as 
as well their particular frequencies (?), and many analyti- 
cal and numerical studies have been presented, as recently 
summarised by (?, hereafter, Paper 1). However, the growth 
of the oscillations will be damped by a magnetic field. A 
transverse field not only reduces the immediate post shock 
temperature and compression but could also completely sta- 
bilise the cooling layer. The first question is: how strong does 
the field have to be to ensure stability? This question has 
been answered for the important case of fast shocks into 
atomic gas (??). Here, we extend these results to include 
power-law cooling functions pertaining to other interstellar 
conditions: a molecular gas and a medium in which the en- 
ergy levels relevant to the cooling are maintained in local 
thermodynamic equilibrium. The conclusions take on added 
significance through the latest campaigns to explore infrared 
and submillimetre regimes with spectroscopic methods (e.g. 
with Spitzer and Herschel). 
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A linear analysis of plane-parallel radiative shocks with 
a transverse field was performed by ?. They assumed a spe- 
cific heat ratio (7) of 5/3 and a cooling function (A) oc p 2 T a , 
where p is the density and T is the temperature. Only the in- 
dex a determines the stability in the absence of a magnetic 
field provided the shock is strong. The latter assumption 
eliminates the Mach number, M, from the problem. Here, 
M = Ui n /c s where u; n is the shock speed and c s is the 
upstream sound speed. 

On the other hand, the Alfven number M a — Ui n /v a , 
where v a is the upstream Alfven speed, is a second variable 
when a significant magnetic field is present. (?, hereafter 
TD93) found, as an example, that even a quite weak field 
(M a < 8) will suppress the growth of all modes of oscillation 
for a > 0. For a = 0.5, even weaker fields (M a < 33) are 
enough to stabilise a shock. Their numerical study agreed 
with the linear analysis and also revealed that for shock 
speeds v s < 160 km s -1 , radiative shocks occurring in in- 
terstellar gas with rin < 0.4 cm -3 may be magnetically sta- 
bilised. Further simulations by ? however, also demonstrated 
that the typical interstellar field may not be sufficient to sta- 
bilise shocks if a < 0. The fundamental mode can generate 
very large amplitude oscillations. Moreover, even when the 
fundamental is stabilised, the overtones can still produce 
higher frequency oscillations of substantial amplitude. 

The general linear analysis for hydrodynamic flows was 
extended in Paper 1 to include molecular shocks. The de- 
pendence on three parameters was considered: a, f3 and 7 
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where j3 is the density dependence of the cooling: 



Aoc / T a . 



(1) 



A strong dependence on 7 was found, with the regime of 
overstability significantly reduced for molecular shocks. In 
particular the fundamental mode grows only for a < —0.24) 
in the molecular equivalent (7 = 7/5 and f3 = 2) of the 
atomic case (for which a < 0.38 is required). However, the 
overtones are more significant in molecular shocks with, for 
example, the first overtone growing for a < 0.66. Building 
on this work, we introduce the Alfven number, here as the 
fourth parameter, to determine the influence of the magnetic 
field. 

In the molecular shocks discussed here, we presume that 
the fraction of mass in ions is sufficiently high so that ion- 
magnetosonic waves cannot propagate ahead of the shock 
front. Otherwise, if the fraction of ions would be low, the 
magnetic field will drift through the neutral molecules, a 
process termed ambipolar diffusion. The friction between 
the ions and neutrals may then replace the molecular viscos- 
ity, changing the nature of the shock front. We discuss this 
restriction in § 01 Other dynamic instabilities are then pos- 
sible in such Continuous Shocks (as well as in Jump Shocks, 
see Paper 1 for a summary). The Jump Shocks studied here 
thus require a source of ionisation that does not cause whole- 
sale molecule dissociation. This may occur near the edges 
of molecular clouds partially exposed to the ultraviolet from 
nearby OB associations, as proposed for the HH 90/91 shock 
(?). The ions can be dust grains, molecules or atoms, accord- 
ing to the physical conditions. 

In this paper, we extend the formalism of TD93. Start- 
ing with the steady state radiative shock structure, we per- 
turb the shock velocity and determine the growth or damp- 
ing rate in addition to the oscillatory period. We study the 
one dimensional case while ignoring thermal conduction as 
well as radiation transfer. In the hydrodynamic 
tionary wall of infinite density is assumed to lie downstream 
of the shock front. In the MHD analysis, the density in 
the downstream gas rises towards a finite constant value. 
In order to make the problem tractable, we not only en- 
sure that the temperature continues to fall, reaching zero 
Kelvin, but that the temperature perturbation also van- 
ishes at the downstream boundary. Then, waves escaping 
downstream are constrained to be pure Alfven waves, lead- 
ing to a straightforward boundary condition. In this respect, 
it should be noted that the published versions of equation 
(A27) of TD93 (where the minus sign should be replaced by 
a +) and equation (2) of ? (where the conditionals should 
be reversed) are incorrect. 

In Paper 1, we undertook a linear stability analysis 
similar to ? while taking a general 7 so as to encompass 
molecular shocks. By also generalising the density index, we 
take into account different physical situations. In particu- 
lar, P = 1 corresponds to cases where the energy levels of 
the atoms or molecules which provide the dominant cool- 
ing are in local thermodynamic equilibrium. For example, 
this can correspond to H2 cooling in warm gas for densities 
above ~ 10 4 cm -3 provided hydrogen atoms act as the main 
collision partner. The specific heat ratio 7 = 5/3 remains ap- 
propriate for dissociative shocks. The value 7 = 7/5 applies 
to a pure H2 gas while, more realistically for the interstellar 
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Figure 1. A sketch of the steady shock configuration. The shock 
jump conditions are defined by equations ll(j| [T71 and [TH1 the to- 
tal compression, —Wf, is given by equation 1221 Note that under 
strong shock conditions, the upstream thermal pressure is equated 
to zero. 



medium, 7 = 10/7 accounts for the addition of ten per cent 
helium atoms. 



2 FORMULATION OF THE PROBLEM 
2.1 The equations 

We consider a reference frame in which an ambient gas of 
uniform density p a and supersonic speed Ui„ passes through 
a stationary shock front at x = x s from upstream i.e. to the 
right in Fig.Q as defined by ?. The gas cools and collapses 
before reaching a cold shell at x = 0, also of uniform density. 
The approach speed m n is defined as a positive quantity. 
Therefore, the pre-shock gas velocity is v = —v,i n . In this 
frame, the material in the shell is moving downstream, thus 
making space for the accumulation of new material. 

Following TD93, we take a perpendicular magnetic 
field, B a , which is frozen into the gas and so increases in 
proportion to the density. The shock Alfven (Mach) number 
is defined as M a = Ui n /v a where v a is the Alfven speed in 
the ambient medium. The magnetic pressure in the flow is 
then 



Pb = bp\ 

where the constant b is 



b = -5 — 2 = v a /2, 



(2) 



(3) 



and we define the useful parameter 9 (the parameter f3 of 
TD93) as 

?2 



B, 



1/(2 Mi), 



(4) 



8np a u'f n 

which is the ratio of the magnetic pressure to the upstream 
ram pressure. 
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Figure 2. A graphical summary of results for the stability. The growth/damping rates, <5 r plotted against a for the fundamental, first 
and second mode (from left to right), the indicated different cases of 7 (from top to bottom), different (3 (thick: /3 = 2; thin; (3=1) and 
three magnetic field strengths, indicated by the line styles shown in the top left-hand panel by the shock Alfven numbers, M a . 



The one dimensional hydrodynamical equations are 
(e.g. TD93) 



dp d(pv) = Q 
dt dx 
' dv dv\ dP 
dt dx ) dx 



and 



de ^ d(ev) _ d(Pv) ^ 



dt dx 



dx 



(5) 
(6) 

(7) 



where A is the radiated energy loss per unit volume and 
P is the total (magnetic+thermal) pressure related to the 
internal energy by 



1 - 



+ 7;P V 



(8) 



7 - 1 ' V 7 - 1 / 8vr 2 
Equations JSJ, & Q refer to the conservation of mass 



momentum and energy, respectively. Eliminating e, 
dP dP r n L a ,„ \ I dv 

-( 7 -l)A/-"(P-6p 2 r, (9) 

where A is a constant and we have used the ideal gas law to 
eliminate temperature. 



2.2 The steady state 

The steady state solution is denoted by the subscript 0. 
Equations JSJ and © are integrated to yield 



PQV0 = —paUin 

Pa = paUi„(vo + Uin) + bpa 



(10) 

(11) 



Equations J I (It & 1111 are substituted in equation © which 
results in 



dx 



C(—p a Ui„) 9 1 [— Vq — UinVQ + p a b 


' "in, _ "0 N 
V v u in ) 
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v [ v + 7(^0 + Uin) + Pab 


yS~ + v l )_ 





, (12) 
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Figure 3. A graphical summary of results for the frequency. The frequency, Si, plotted here against a has been normalised by the factor 
l/(2n + 1) where n is the mode number (fundamental being n = 0), which demonstrates that the modes resemble that of a half-open 
organ pipe. Displayed are results for the fundamental, first and second mode (from left to right), the indicated different cases of 7 (from 
top to bottom), different f3 (thick: /3 = 2; thin; /3 = 1) and three magnetic field strengths, indicated by the line styles shown in the top 
left-hand panel by the shock Alfven numbers, M a . 



where 

C= (7-1)* A 

This equation can be integrated to determine the velocity 
through the shocked layer. We now introduce the following 
variables: 



c = - 



"ill 



(13) 
(14) 



Equations and O lead to 

-{-wf \w + 7(1 + w + 6) + i^jli 



(15) 



dw " v?°T*Cpl- x x a {-w - - dw + |) a ' 
The jump condition across the steady shock front can be 



written as (? 

Po(xs) = 
v (x 3 ) = 

Po(x s ) 



Paq, 

Ui 



p a uf n ( 1 - - + I 

q 



(16) 
(17) 

(18) 



where the shock front compression is 

• - m>- 

U = -(7-1 + 27^), 

V = [^470(70 - (7 - 1)) + (7 - I) 2 + 160"] . 

The boundary condition at x = is determined from the 
fact that the temperature is zero: the total pressure is equal 
to the magnetic pressure and the density remains finite: 

PoOe = 0) = P B (x = 0). (20) 

If we denote the velocity of the shell as vq(x = 0), and also 
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Figure 5. The growth rate for various overtones in the presence 
of a strong transverse magnetic field for the cases a = 0.3 (upper 
panel) and a = —0.5 (lower panel) on taking /3 = 2, 7 = 5/3 and 
M a = 1.5. 



Figure 4. The loci of neutral stability for the indicated values of 
7, modes and = 2 represented by thick lines and j3 = 1 by thin 
lines, log (M a ) is plotted against a). The region of instability is 
above the lines in all cases. 



define w f = we can then write from ||5J and HH 

the expression 



PaUi„(l + w f + 6) 



PaUi 



The physical solution to the above cubic equation is 



Wf 



^/(6 2 +46) 



Note that —1/wf is the total compression, as indicated in 
Fig- El 



2.3 The set of linear equations 

The shock wave is now perturbed by 

or = v ° ie ' (23) 

where a = tjR+iai is the frequency and v„i is a real quantity. 
The position of the shock may be represented as the real part 



of 

x a = x 30 + x sl e' Tt (24) 
where a; s i = v al /a. Considering only the terms up to first 



(25) 
(26) 
(27) 





order: 




(21) 


X s 


X f Xsl 
X B V X a o 




2i 


XsO V X s0 




dx 


(22) 


d(L 
dt 


xx s -i_ae at 

1 ,s() 



p 
p 



p (0 + Pi(0e°' 



(28) 
(29) 
(30) 



All the quantities with subscript 1 represent the small per- 
turbed factors. 

When the shock is moving, we consider the frame of 
reference in which the shock is stationary, i.e., the upstream 
velocity is now — m n — t^ie 17 '. The boundary conditions at 
the moving shock wave are obtained by taking the deriva- 
tives of equations 116H . I17H and 118H with respect to Ui n to 
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obtain 
pi(x = Xs) 

P\(X = Xs) 
Vi(x = Xs) 



2 Pa d(l/q) 



dO 



2v a \p a U in yi - 



} d(l/q) \ 

de ) 



vsi 29 



de 



+ V s 



(31) 
(32) 
(33) 



where the last term in 133H arises because ol the transfor- 
mation of coordinates to the shock frame. 

We will also need the derivative of q which is, from 
equation 1191 . 



d(l/q) 



dO 



1 + 



4 + 7(2761- (7-I)) 



2-1 



.2-7(7-1). 
We then transform to the following variables: 
_ x s oapi 

Vslpa 

Pi 



VslPaUin 
Vl 

Xsoa 



(34) 

(35) 
(36) 
(37) 
(38) 



At the shock front, from 13111 - 1331 1. the new variables take 
the following values: 



w s ) 


= 1, 


w s ) 


= 256q 


W s ) 


-,(, 


Ws) 


- (V 



a d(l/g) 

de '■ 



+ 



de J ' 



d (l/q) 

de 



+ 1 



(39) 
(40) 

(41) 
(42) 



Substituting gg), and into ©, © and ©> the 
fluid equations can be rewritten as 



. df; w dC, 



1 d V C 



' dw 



8 dw w 2 w dw 

„2 



d£ dr\ w dir 

-? + — V + w— +77 j-C - w — 

dw aw dw 



where 
D = 

E = 



D + E 



= 0, 


(43) 


= 0, 


(44) 


= F 


(45) 



r df dir w 

-§ + 0— 7T + w- \-rj - 

dw dw 



+ 



d£ w d( 
dw 5 dw 

(2-7)6 



+ 7T7 



010 



i» + 7(l + ^ + g)+ { - 2 ~^ 6 



2(2 - 7 )C6» 



cwr - 08 - q)(1 + 10 + 6>)^ + (/3 + a)|i 
ow 



w + 1 + e — - 



The quantities n and r\ are complex eigenfunctions. We 



employ the subscript r to denote the real component and i 
stands for the imaginary part for each of the above quanti- 
ties. The quantity 5 is a complex number with the sign of the 
real part, S r , indicating the instability (+ve value) or stabil- 
ity (-ve value) of a mode. The quantity 5, is interpreted as 
the eigenfrequency (in units of («i„/i s o)). The equations for 
the real and complex quantities are provided in AppendixIXl 



2.4 The integration process 

While integrating the above equations we find that there 
are numerical problems when a < 0.2 as the temperature 
approaches zero. To solve this, we follow a prescription sug- 
gested by TD93 which is to introduce a break in the power 
law cooling at a temperature T c , at a velocity w c , leaving a 
small region near the shell where the cooling function will 
have a = 0.5. Therefore, the cooling function is split into 
two regimes such that 



A = V" (P ~ bpT 



A = A c p< 



{3—a. c 



(p - bp 2 



at w < w c 
at w > w c 



(46) 
(47) 



Here, A c = A(P C - bp c ) a ~ a " and we choose A c = 0.001 to 
be consistent with the results of TD93. To obtain the shock 
length, we integrate equation 1151 from £ = 1 to £ = and 
take into account the two component cooling function. 

The boundary condition at the shell for a moving shock 
under the influence of a transverse magnetic field is derived 
in Appendix^] It is more complicated than in Paper I where 
the gas settles on the wall. In the present scenario, the ve- 
locity of the gas does not vanish since the magnetic pressure 
limits the compression. The boundary condition is 



'20 



+ n\ = 0, 



(48) 



similar to TD93 but correcting a typographical error in their 
equation (A27). 

For the various combinations of the four free parame- 
ters, a, P, 7 and e, the quantities S r and Si are determined 
by imposing the boundary condition at the shell boundary. 
We solved the differential equations employing a fourth or- 
der Runge-Kutta technique for trial values of S r and <5, se- 
lected from a grid of points uniformly covering the complex 
plane. The combinations that come closest to satisfying the 
boundary condition for each mode determine a new set of 
grid points with a higher resolution. 



3 GENERAL RESULTS 

We have calculated the growth rates (S r ) and eigenfrequen- 
cies (Si) for the various cases of a, (3 and 7 with weak 
(M a = 30), intermediate (M a = 10) and strong (M a = 3) 
magnetic fields for the fundamental mode as well as the first 
three overtones. These results are presented in Tables IClT - 
IC8I Equivalent tables for the case 7 = 10/7 have been con- 
structed but are not presented here. We plot the growth 
rates against a in Fig. [5] for various modes. The major re- 
sults are as follows. 

• We have accurately reproduced the results of TD93 for 
the specific case (3 = 2 and 7 = 5/3 ( Tables IC2lfc Table IC6|l 
as a specific case in our parametric space. 
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• The fact that the instability regime becomes increas- 
ingly restricted with an increase in the transverse magnetic 
field holds for all values of a, (3 and 7. 

• The fundamental remains the most restricted mode 
to be unstable. The field has a similar stabilising effect on 
all the modes. 

• The lower value for the density index /3 provides more 
stability. However, for high magnetic fields (solid lines in 
Fig. |2J, the stability criterion becomes independent of (3. 
This is interpreted as the cushion effect of the field which 
limits the compression. 

• With decreasing 7, the instability regime is more re- 
stricted for all magnetic field strengths. Whether the molec- 
ular gas is diatomic (7 = 7/5) or contains helium (7 = 10/7) 
does not significantly alter the instability regime. 

• The frequency of the fundamental decreases strongly 
as a decreases for all magnetic field strengths. The frequency 
of the overtones generally tend to a constant value as a 
decreases. Note that in Fig. [3] the frequencies of the modes 
have been normalised by a factor of 2n + 1 where n is the 
mode number for display purposes and so demonstrates that, 
for all magnetic field strengths, the resonant frequencies are 
analogous to the acoustical modes in a half-open organ pipe. 

The loci of neutral stability in the log(M a ) - a space, 
which separate stable and unstable regimes (for fixed /3 and 
7) , are displayed in Fig. [I] Whereas an Alfven number of 
8 will suppress the growth of all modes of oscillation for 
a > in the atomic shocks (7 = 5/3, /3 — 2), we require just 
M a < 20 for the equivalent molecular case. Hence molecular 
shocks tend to be considerably more stable. On the other 
hand, the Alfven number of 8 would be sufficient to stabilise 
even the first overtone for molecular shocks with a = —0.5. 

Given the trends in Fig. 2] it is not clear that higher 
order modes possess wider instability ranges. In fact, from 
inspection of the tables (in the appendix), the growth rate 
tends to be lower for higher overtones for a > 0.0 especially 
as the magnetic field strength increases, whereas for a < 0, 
however, there is an increasing trend in the growth rate. 
Here, we have investigated these trends for low M a shocks to 
determine if the magnetic field will stabilise shocks entirely. 
As shown in Fig. the growth rate continues to decrease as 
the mode number is increased for a = 0.3 and M a — 1.5, tak- 
ing an atomic shock for illustration. In contrast, the growth 
rate displays the opposite trend for the case with negative 
q. In other words, the magnetic field may not provide com- 
plete stability for cooling with negative a. However, the sig- 
nificance of high order modes are probably highly restricted 
and a multi-dimensional analysis is necessary to determine 
any observational consequences. 

The present analysis applies to both atomic and molec- 
ular shocks and assumes that the magnetic field is parallel to 
the shock front. However, in a sufficiently strong shock the 
molecules will dissociate immediately following the shock 
jump. In fact, dissociative cooling may dominate the cool- 
ing function. Given a molecular medium with Alfven speed 
of 1 km s" 1 (?) and a speed limit for molecular shocks 
of 40 km s _1 at low densities (under ~ 10 4 cm -3 ) and 
~ 24 km s _1 at high cloud densities (?), the maximum 
Alfven number for molecular shocks with 7 = 7/5 is then in 
the range 24 - 40. 

A maximum value for the magnetic field should also 
exist, above which ambipolar diffusion dominates the shock 



physics. The limit is more difficult to calculate since it de- 
pends on the ionisation fraction, density and cooling func- 
tion. In fact, the limitation is much stronger on the ion 
fraction rather than the magnetic field. In molecular clouds 
which are well shielded from the extreme ultraviolet, the 
ion fraction, Xi is extremely low and ambipolar diffusion 
dominates. On the other hand, the ion fraction may reach 
high values near massive stars, stellar ouflows or within bow 
shocks with high-speed apices, supplying ionising radiation. 
As estimated by ?, ion fractions exceeding ~ 10~ 5 B-3/ng 
(H 2 cooling) or ~ 2 10" 6 B- 3 /n\ /2 (H 2 cooling) are re- 
quired to ensure a J-type shock where B_ 3 = B/(10~ 3 G) 
and n% = n/(10 6 cm -3 ). Thus, the Alfven speed must not 
exceed ~ 5 W 5 neX (H2O cooling) km s _1 (H2O cooling) or 
~ 10 5 x km s" 1 (H 2 cooling). 



4 CONCLUSIONS: MOLECULAR SHOCK 

We have found the stable and unstable regions for radiative 
shocks in the a, (3 and 7 parameter space for various mag- 
netic field strengths and in the field-Q parameter space for 
specific [3 and 7 values. 

For the fundamental mode, molecular shocks are con- 
siderably more stable than atomic shocks. For the overtones, 
molecular shocks are only moderately more stable. The mag- 
netic field, however, has a strong stabilising influence on all 
these modes, decreasing all the growth rates by similar fac- 
tors (see Fig. |2J. 

The field strength enters the stability condition through 
the shock Alfven number. All cooling functions in which the 
cooling increases with temperature are stable for M a < 20. 
However, for M a > 22, the first overtone can be unsta- 
ble, potentially leading to oscillations with non-linear am- 
plitudes. Hence, we need to estimate the Alfven number for 
typical shocks in various environments. 

Unfortunately, in molecular clouds, the Alfven speed is 
a very uncertain parameter. It is even difficult to get an 
objective view on the magnetic field strength alone (?) al- 
though derived Alfven speeds generally lie within the range 
0.5-5 km s _1 (?). Jump shocks followed by cooling zones 
can arise without molecular dissociation for shock speeds as 
high as 40 km s _1 at low densities but ~ 24 km s _1 at high 
cloud densities (?). As discussed in § a minimum ion frac- 
tion is also required to maintain a frozen-in magnetic field. 
Hence, a quite fast shock in a diffuse molecular cloud is the 
most favourable state for the appearance of oscillations. 

The temperature dependence of molecular cooling func- 
tions generally correspond to values of a well above zero and 
so indicates stable shocks. However, besides temperature 
and density, shock cooling may depend on the chemistry. 
The rate of formation of trace molecules within the cooling 
layer can increase the cooling rate as the temperature falls. 
In particular, the propagation of warming but non-ionising 
shocks into cool atomic gas may show trace molecule forma- 
tion and increasing cooling within a gas with 7 = 5/3. How- 
ever, because of the low shock speeds, these shocks would 
be stabilised if the magnetic field were transverse. Hence, it 
can be that oblique non-ionising shock waves into atomic gas 
represent the rather restrictive conditions for shock instabil- 
ity within the cool components of the interstellar medium. 

Although it is well known how an oblique field alters 
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the nature of a steady shock, it is not clear how it alters d( r _ £5r _ d£_ CrSr 5r_ dr/ r _ rj r 5 r 

the stability conditions. It should also be noted that the dw w 3 dw w w 2 dw w 3 

magnetic field may help stabilise shocks even when the field _ S r (Cr5 r + Ci&) _j_ d£_ _ df]i 

is parallel to the shock provided the Alfven number is less wS 2 dw w w 2 dw 

than the Mach number, a condition expected to be satisfied ViSj $j{CiS r — Crj») (A5) 

in molecular clouds (?). Such an analysis, as well as multi- w 3 w8 2 

dimensional numerical simulations of radiative shocks, are d£i_ _ C^i _ d^ C,i8 r S r dr\j ^ r]iS r 

still to be performed. dw w 3 dw w w 2 dw w 3 

wS 2 dw w w 2 dw 

w 3 w5 2 
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APPENDIX B: BOUNDARY CONDITION AT 
THE SHELL 

Following TD93, we exploit the condition that the temper- 
ature tends to zero at the shell. This implies that the only 
pressure at the shell is due to the magnetic field. To extract 
the boundary condition, we represent the physical quantities 

as 



p = po(x) + pi(x)e u T ' 
P = Po(x) + Pi(x)e l(/ * +fea;) 
v = vo(x) + vi(x)e l ^ t+kx ^ 



(Bl) 
(B2) 
(B3) 



APPENDIX A: 
EQUATIONS 



THE DIFFERENTIAL 



^,From Eqs. Ij36^ . I|37p and l|38|l . we get six coupled first order 
equations which are 



dr/r 
dw 



dr) t 
dw 



dn r 
dw 



+ 



2^5 2 - 2r) r 8 2 



(l + w + 9-^)S 2 

JTTrS 2 + (W 2 + Hi2_!^)( Cr(5r + Q iS . 



( w + j( w + i + e) + ^yii)s 2 



dw 



(w + j(w + 1 + 1 



, (2-7)9 > 



5 r 
+ — 

S 2 



where k is the wavenumber, which represents the direction of 
the outgoing wave towards the shell and / is the frequency. 
The above equations are substituted in © and © to yield 



rS 2 + (S r ( r + (^±Hl« _ w{/3 - a )(l + w + 9)) { = -[ 



I 

k 

I 

k 



f povivo + Pi 



povi 

POVl + PIVQ 

pi 



Equating the above expressions leads to 



Pi = 



2 2 
Pi 



(Al) 



The magnetic pressure is given by 



OLTYiS 2 + (8 r (i — Si^r) 



— w(f3 — a)(l + w + i 



8tt 



(B4) 
(B5) 

(B6) 

(B7) 



+ ■ 



(l + w + e-^)5 2 
1 ^8 2 + (w 2 + 2 J 2 ^)(08 r 



CrSi) 



We now express the magnetic field in terms of density for 
example, (?) as 



(w + j(w + 1 + 9) + 



(2-7)9 



dw 



5i(r) r + 7Tr) + S r (rji + 7T») 

(w + 7 («, + 1 + 9) + ^yR) 



B_ 



£_ 

Pa 



(B8) 



5i_ 

5 2 



(A2)where B and p are the post-shock magnetic field and density. 
Substituting <B2> in 1B7I I and defining the Alfven velocity 

as v A = J £7, we get 



-£ (SrTjr - 5jJ]i) rig TJr_ dr/r 

w w dw w dw 

W)(Cr<5r + (iSi) 



Pi 



vaPi 



(B9) 



S 2 



(A3) 



Equations IIB6II and lB9t yield the boundary condition at 
the shell as, 



d-Kj _ {5jT} r + S r T]i) d£ Tfr ^ (k]^ 

dw w dw w dw 



w 

w(£i8 r 



CrSi) 



d 2 



(A4) 



Pi 



-v avi pa 



(BIO) 



The minus sign indicates the direction of the gas flow. Note 
that the perturbed quantities are function of x, which itself 
is a function of £ from 125H . The expression 1B101 can be 
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written in terms of the non-dimensional quantities as 




(Bll) 



The above quantities are evaluated at the shell boundary or 
equivalently at w = w/. 
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Table CI. Growth rates (<S r ) for several modes and magnetic field strengths, as measured by the Alfven number M a , for 7 = 



5/3, /3 = 1. 







fundamental 






first overtone 




a 


M a = 3 


M a = 10 


M a = 30 


Ma = 3 


M a = 10 


Ma = 30 


-1 


- 0.02942 


0.04903 


0.08965 


0.06387 


0.11420 


0.14986 


-0.5 


- 0.08446 


0.00567 


0.03646 


- 0.03742 


0.03890 


0.07482 


0.0 


- 0.18538 


- 0.05842 


- 0.02296 


- 0.20521 


- 0.06762 


- 0.01974 


0.3 


- 0.29720 


- 0.11769 


- 0.07122 


- 0.38342 


- 0.16282 


- 0.09245 


0.5 


- 0.42465 


- 0.17862 


- 0.11462 


- 0.57743 


- 0.25705 


- 0.15414 







second overtone 




third overtone 




a 


M a = 3 


M a = 10 


M a = 30 


M a = 3 


M a = 10 


M a = 30 


-1 


0.11050 


0.13937 


0.16695 


0.13817 


0.16075 


0.19142 


-0.5 


- 0.00862 


0.04990 


0.07818 


0.00809 


0.06541 


0.09936 


0.0 


- 0.20728 


- 0.07856 


- 0.03713 


- 0.21087 


- 0.07366 


- 0.02433 


0.3 


- 0.42137 


- 0.19620 


- 0.12957 


- 0.45058 


- 0.20419 


- 0.12893 


0.5 


- 0.65377 


- 0.31342 


- 0.21165 


- 0.70913 


- 0.33717 


- 0.22760 



Table C2. Growth rates (<5 r ) for several modes and magnetic field strengths for 7 = 5/3, (3 = 2. This Table can be compared to Table 1 
of TD93. 







fundamental 






first overtone 




a 


M a = 3 


M a = 10 


M a = 30 


M a = 3 


M a = 10 


M a = 30 


-1 


- 0.03549 


0.04610 


0.09351 


0.06696 


0.13981 


0.19781 


-0.5 


- 0.08877 


0.01342 


0.04383 


- 0.02921 


0.08413 


0.13291 


0.0 


- 0.17983 


- 0.02853 


0.00702 


- 0.18560 


0.00496 


0.06983 


0.3 


- 0.27889 


- 0.06263 


- 0.01806 


- 0.35272 


- 0.06229 


0.02479 


0.5 


- 0.39054 


- 0.09283 


- 0.03697 


- 0.53889 


- 0.12623 


- 0.00998 







second overtone 




third overtone 




a 


M a = 3 


M a = 10 


M a = 30 


M a = 3 


M a = 10 


M a = 30 


-1 


0.11374 


0.17466 


0.22738 


0.14253 


0.19826 


0.25109 


-0.5 


0.00067 


0.10719 


0.16104 


0.01897 


0.12425 


0.18399 


0.0 


- 0.18654 


0.00890 


0.08662 


- 0.18846 


0.01377 


0.10306 


0.3 


- 0.38992 


- 0.08046 


0.03131 


- 0.41714 


- 0.09122 


0.04081 


0.5 


- 0.61605 


- 0.17289 


- 0.01352 


- 0.67020 


- 0.20545 


- 0.01179 



Table C3. Growth rates (5 r ) for several modes and magnetic field strengths for 7 = 7/5, (3 = 1. 







fundamental 






first overtone 




a 


Ma = 3 


Ma = 10 


M a = 30 


Ma = 3 


Ma = 10 


Ma = 30 


-1 


- 0.12574 


- 0.03009 


0.00029 


- 0.04055 


0.02876 


0.06090 


-0.5 


- 0.18312 


- 0.06286 


- 0.02743 


- 0.13786 


- 0.02838 


0.01337 


0.0 


- 0.28865 


- 0.11479 


- 0.06609 


- 0.30446 


- 0.11369 


- 0.05206 


0.3 


- 0.40691 


- 0.16801 


- 0.10022 


- 0.48460 


- 0.19620 


- 0.10778 


0.5 


- 0.53663 


- 0.22589 


- 0.13324 


- 0.67452 


- 0.28096 


- 0.16065 








second overtone 




third overtone 




a 


M a = 3 


M a = 10 


M a = 30 


Ma = 3 


M a = 10 


M a = 30 


-1 


0.00362 


0.04071 


0.05939 


0.02880 


0.06460 


0.08899 


-0.5 


- 0.11048 


- 0.02764 


0.00114 


- 0.09556 


- 0.09275 


0.02439 


0.0 


- 0.30800 


- 0.13002 


- 0.07824 


- 0.31358 


- 0.12335 


- 0.06976 


0.3 


- 0.52493 


- 0.22824 


- 0.14361 


- 0.55789 


- 0.23731 


- 0.15350 


0.5 


- 0.74988 


- 0.32688 


- 0.20005 


- 0.80390 


- 0.35412 


- 0.22737 
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Table C4. Growth rates (S r ) for several modes and magnetic field strengths for 7 = 7/5, (3 = 2. 







fundamental 






first overtone 




a 


Ma = 3 


M a = 10 


Ma = 30 


Ma = 3 


Ma = 10 


Ma = 30 


-1 


- 0.f2850 


- 0.02761 


0.01627 


- 0.03781 


0.04672 


0.08856 


-0.5 


- 0.18572 


- 0.04914 


- 0.01205 


- 0.13274 


0.00162 


0.05346 


0.0 


- 0.28557 


- 0.08540 


- 0.03462 


- 0.29240 


- 0.06666 


0.01054 


0.3 


- 0.39682 


- 0.11950 


- 0.05177 


- 0.46628 


- 0.13231 


- 0.02255 


0.5 


- 0.51900 


- 0.15418 


- 0.06567 


- 0.65304 


- 0.20231 


- 0.05032 







second over 


tone 




third overtone 




a 


M a = 3 


Ma = 10 


M a = 30 


M a = 3 


M a = 10 


M a = 30 


-1 


0.00628 


0.06857 


0.10574 


0.03269 


0.09028 


0.13057 


-0.5 


- 0.10479 


0.01285 


0.06522 


- 0.08846 


0.02882 


0.08690 


0.0 


- 0.29458 


- 0.07317 


0.01163 


- 0.29851 


- 0.06802 


0.02653 


0.3 


- 0.50449 


- 0.15932 


- 0.03224 


- 0.53521 


- 0.16798 


- 0.02552 


0.5 


- 0.72683 


- 0.25388 


- 0.07193 


- 0.77992 


- 0.27967 


- 0.07602 



Table C5. Eigenfrequencies (Si) for several modes and magnetic field strengths for 7 = 5/3 and (3=1. 







fundamental 






first overtone 




a 


M a = 3 


M a = 10 


M a = 30 


Ma = 3 


M a = 10 


M a = 30 


-1 


0.44847 


0.33935 


0.30432 


1.45520 


1.03328 


0.92059 


-0.5 


0.50430 


0.37630 


0.34126 


1.50546 


1.03500 


0.91030 


0.0 


0.58627 


0.41387 


0.36464 


1.59283 


1.03047 


0.87269 


0.3 


0.67073 


0.44723 


0.37841 


1.68773 


1.02950 


0.83627 


0.5 


0.75300 


0.48020 


0.38821 


1.76664 


1.02936 


0.80303 



second overtone third overtone 



a 


Ma = 3 


M a = 10 


M a = 30 


Ma = 3 


M a = 10 


M a = 30 


-1 


2.50463 


1.77852 


1.59824 


3.55657 


2.51123 


2.25140 


-0.5 


2.54768 


1.74493 


1.53713 


3.59392 


2.44672 


2.14906 


0.0 


2.63290 


1.69280 


1.43223 


3.67699 


2.35191 


1.98756 


0.3 


2.72930 


1.65092 


1.33791 


3.77286 


2.27393 


1.84544 


0.5 


2.79538 


1.60990 


1.24986 


3.82144 


2.19758 


1.71171 



Table C6. Eigenfrequencies (Si) for several modes and magnetic field strengths for 7 = 5/3 and (3 = 2. 







fundamental 






first overtone 




a 


M a = 3 


M a = 10 


M a = 30 


M a = 3 


M a = 10 


Ma = 30 


-1 


0.43031 


0.27658 


0.25635 


1.45075 


1.03787 


0.94698 


-0.5 


0.47773 


0.32531 


0.28242 


1.49626 


1.05480 


0.95981 


0.0 


0.55034 


0.35176 


0.30317 


1.58012 


1.06036 


0.94929 


0.3 


0.62587 


0.37345 


0.31234 


1.67378 


1.07176 


0.93804 


0.5 


0.70256 


0.39530 


0.31871 


1.75753 


1.09095 


0.92967 








second overtone 




third overtone 




a 


M a = 3 


M a = 10 


Ma = 30 


M a = 3 


M a = 10 


M a = 30 


-1 


2.50320 


1.81950 


1.68275 


3.55754 


2.59039 


2.40284 


-0.5 


2.54284 


1.80897 


1.66530 


3.59230 


2.55669 


2.36004 


0.0 


2.62573 


1.78459 


1.62070 


3.67348 


2.50237 


2.28419 


0.3 


2.72132 


1.77329 


1.58450 


3.76880 


2.46542 


2.22368 


0.5 


2.79310 


1.77345 


1.55688 


3.82473 


2.43991 


2.17669 
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Table C7. Eigenfrequcncies (Si) for several modes and varying magnetic field strengths for 7 = 7/5 and f3 = 1. 







fundamental 






first overtone 




a 


M a = 3 


M a = 10 


M a = 30 


M a = 3 


M a = 10 


M a = 30 


-1 


0.43865 


0.27943 


0.23392 


1.41341 


0.84427 


0.69527 


-0.5 


0.48338 


0.30069 


0.25072 


1.46172 


0.84617 


0.68237 


0.0 


0.55590 


0.32986 


0.26344 


1.55387 


0.85161 


0.65657 


0.3 


0.62987 


0.35788 


0.27225 


1.64946 


0.86141 


0.63397 


0.5 


0.69568 


0.38442 


0.27884 


1.71604 


0.86861 


0.61193 








second overtone 




third overtone 




a 


M a = 3 


M a = 10 


M a = 30 


M a = 3 


M a = 10 


M a = 30 


-1 


2.43144 


1.44978 


1.19709 


3.45026 


2.04889 


1.69333 


-0.5 


2.47796 


1.42644 


1.14496 


3.49517 


2.00664 


1.61494 


0.0 


2.57765 


1.39673 


1.06321 


3.60095 


1.95155 


1.49445 


0.3 


2.68125 


1.37676 


0.99033 


3.70963 


1.90977 


1.38357 


0.5 


2.73264 


1.35376 


0.92053 


3.74108 


1.85977 


1.27176 



Table C8. Eigenfrequcncies (Si) for the various modes of parameters 7 = 7/5 and /3 = 2 for various magnetic fields. 







fundamental 






first overtone 




a 


M a = 3 


M a = 10 


M a = 30 


M a = 3 


M a = 10 


M a = 30 


-1 


0.42613 


0.24811 


0.19424 


1.41061 


0.84868 


0.71088 


-0.5 


0.46608 


0.26937 


0.21581 


1.45638 


0.85532 


0.71608 


0.0 


0.53258 


0.29198 


0.22668 


1.54505 


0.86521 


0.70578 


0.3 


0.60136 


0.31426 


0.23328 


1.63895 


0.88205 


0.69826 


0.5 


0.66468 


0.33806 


0.23891 


1.70814 


0.90367 


0.69442 








second overtone 




third overtone 




a 


M a = 3 


M a = 10 


M a = 30 


M a = 3 


M a = 10 


M a = 30 


-1 


2.42930 


1.47335 


1.25896 


3.44982 


2.09150 


1.79406 


-0.5 


2.47396 


1.46111 


1.23842 


3.49352 


2.06317 


1.75504 


0.0 


2.57015 


1.44656 


1.20060 


3.59567 


2.02658 


1.69333 


0.3 


2.67207 


1.44242 


1.17200 


3.70268 


2.00303 


1.64574 


0.5 


2.72750 


1.44027 


1.15139 


3.73955 


1.97825 


1.60978 



